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SIGNAL ENHANCEMENT MODULE 

CROSS-REFERENCE TO RELATED APPLICATION 
[0001] This application claims benefit of the priority of 
U.S. Provisional Application Serial Number 60/448,200 filed 
February 18, 2003 and entitled "Signal Enhancement Module." 

TECHNICAL FIELD 
[0002] The disclosed system relates to electromagnetic 
radiation imaging systems, and more specifically to filtering 
signals produced from an electromagnetic radiation imaging 
system, based on filtering criteria that imply that the 
signals detected do not originate from or pass through the 
object being imaged. 

BACKGROUND 

[0003] Detector systems, such as gamma ray detectors, may be 
sensitive to non-physical events that cannot be discriminated 
on the basis of the signals from a single pixel. These events 
cause noise in the system. 

[0004] The signal generated by an individual gamma ray is 
very small - of the order of a few hundred to a few thousand 
electrons. Therefore, relatively small extraneous signals may 
be detected as gamma rays and erroneously registered. 
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[0005] Vibration in solid-state gamma cameras is one example 
of a phenomenon that may cause the generation of 
electromagnetic noise signals that are detected as one or a 
plurality of events of interest. Analogously, electrical 
discharges, including discharges driven by potential 
differences below 5 volts, may generate non-physical events. 
Impingement of extraneous high-energy particles such as cosmic 
rays may also cause errors. As a result of Compton scattering 
into multiple detection elements, a high-energy gamma ray 
striking the radiation-sensitive surface of a SPERIS can 
directly produce signals in multiple channels. The gamma ray 
can also produce noise signals in neighboring channels through 
coupling by electrical parasitics. 



SUMMARY 

[0006] A technique of filtering events from an 
electromagnetic detector is described. According to one 
aspect, thresholds are set based on expected amounts of 
electromagnetic radiation impinging on a detector. The output 
signal from the detector is filtered based on the threshold, 
to reject signals that are based on undesired stimuli. 
[0007] In one aspect, different sources of noise are 
characterized according to the quantities of data that will be 
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generated by the detector system in response to signals caused 
by each source of noise. 

[0008] In an embodiment, the detector is a segmented detector 
gamma camera, used in the medical imaging application the 
filter is used to determine The present disclosure describes 
filtering from the events detected by a segmented detector, 
pulse counting electromagnetic radiation imaging system or 
SPERIS, events that cause a plurality of individual events 
concurrent in time. The techniques may be especially useful 
for events that cannot be distinguished on the basis of 
primary particle energy, precise time coincidence or the 
geometric relationship between the individual events caused by 
the event . 

[0009] The filter, or signal enhancement module, 
discriminates against these events by rejecting signals that 
exceed definable thresholds in terms of their concurrence in 
time and one or a plurality of the quantity, global density, 
local density and energies of the individual events caused by 
the events. 

DESCRIPTION OF DRAWINGS 
[0010] These and other features and advantages of the 
invention will become more apparent upon reading the following 
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detailed description and upon reference to the accompanying 
drawings . 

[0011] Fig. 1 depicts a decision tree for processing of data 
in the raw data buffer and copy buffer by the raw multiplicity 
filter and the global density filter, including the assignment 
of data to the candidate buffer. 

[0012] Fig. 2 illustrates data processing performed by the 
local density filter. 

[0013] Fig. 3 depicts the decision tree for processing of 
data by the mode/buffer size (MBS) filter. 

[0014] Fig. 4 shows a decision tree of the blink multiplicity 
filter, cosmic multiplicity filter and the filter state 
filter . 

[0015] Fig. 5 illustrates a decision tree of the average 
energy filter. 

[0016] Fig. 6 shows an exemplary electromagnetic radiation 
imaging system. 

DETAILED DESCRIPTION OF THE INVENTION 
[0017] The scope of the techniques disclosed herein is 
intended for use specifically in electromagnetic radiation 
imaging systems, "ERIS", that produce images by detection of 
photons such as X-rays and gamma rays. In addition, these 
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techniques are intended for use with an ERIS that comprise 
segmented detectors and operate in a pulse counting mode. 
Examples of these systems include commercial gamma ray imaging 
systems manufactured by Digirad Corporation for gamma ray 
imaging in nuclear medicine, and some commercial X-ray imaging 
systems . 

[0018] A segmented detector, jsulse counting electromagnetic 
radiation imaging system, or SPERIS 600, is shown in Fig. 6. 
The device typically has four basic subsystems: (1) one or a 
plurality of detector modules 606 that capture the raw signals 
610; the totality of the individual detector modules being 
denoted the detector 605, (2) a signal processing subsystem 
615 that receives information from the detector 605, and 
converts the information from the detector modules into a form 
620 suitable for further processing in a computer, (3) a 
computer 630 that receives the data 62 0 from the signal 
processing subsystem 615, then generates and may enhance one 
or a plurality of final images and may also calculate the 
values of attributes of the object imaged, and (4) one or a 
plurality of output devices 640 such as computer monitors and 
printers that present the final images and other data in 
desired formats. 

[0019] The minimum requirements for pulse counting, segmented 
detector imaging systems are that the system be able (1) to 
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detect the impingement of particles of radiation of the type 
being imaged on one or a plurality of radiation-sensitive 
regions, (2) to determine the position of impingement of one 
or a plurality of particles on said radiation-sensitive 
surface or surfaces. In many cases, it is also desirable to 
determine the time, with a selected precision, at which said 
particle or particles impinge on said radiation- sensitive 
surface or surfaces. 

[0020] These devices are intended to detect a specified kind 
of radiation, e.g., in the disclosed embodiment, gamma rays. 
Typically the detector module 605 of an electromagnetic 
radiation imaging system will also, however, receive energy, 
e.g., particles of radiation with properties other than those 
of the particles whose image is being generated. Signals 
generated by this undesired radiation may appear as noise in 
an image. 

[0021] The detection of any phenomenon that causes the 
imaging system to generate a signal is referred to as an 
event. An event actually caused by the particle or plurality 
of particles of the type and properties being imaged, and 
originating from or passing through the object being imaged, 
is denoted an event of interest "EOI" . Any event other than 
an EOI is considered to be a non-physical event or NPE. 
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[0022] Typically, SPERIS discriminate between individual 
particles on the basis of energy. This function is often 
performed within the detector module by such operations as 
pulse shaping, peak detection and comparison of peak value to 
a threshold energy reference. Some discriminate between a 
plurality of particles on the basis of time coincidence, and 
some discriminate on the basis of the geometric relationship 
between a pair of particles. A particular ERIS may employ any 
combination of these three discrimination criteria. For 
example, in Anger- type gamma ray imaging systems used in 
nuclear medicine, light generated in a continuous single 
crystal sheet scintillator is collected by a plurality of 
photomultiplier tubes. Because Anger cameras employ a single 
sheet scintillator, they are not SPERIS. A gamma ray 
impinging on the scintillator sheet generates a large number 
of photons that spread throughout the crystal, causing all the 
phototubes to detect a signal simultaneously. Signals that do 
not cause simultaneous detection of light by all the 
photomultiplier tubes are treated as non-physical events. 
Similarly, in positron emission tomography or PET, the EOI is 
the simultaneous emission within the patient of two particles 
of identical, well-defined energy and traveling in 
diametrically opposite directions. Phenomena not meeting 
these criteria are rejected as NPE. 
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[0023] Detection modules for SPERIS comprise a plurality of 
individual detector elements, often called pixels or channels. 
Each individual channel produces signals in response to the 
impingement of a particle or particle. For example, the 
detector module of a solid-state, segmented-detector gamma ray 
imaging system used in nuclear medicine may comprise an array 
of solid-state detector elements and may perform some signal 
processing functions. The array typically comprises a large 
number of pixels and may comprise a plurality of smaller 
arrays. The individual detector elements may comprise a 
direct -conversion detector, such as cadmium zinc telluride, 
which generates electrical signals in response to absorption 
of a particle within the material, or an indirect -conversion 
detector such as a scintillator and a photodiode. Note that 
the Anger camera discussed below comprises a plurality of 
photomultiplier tubes, but only one scintillator crystal and 
is therefore not a segmented detector system 

[0024] The disclosed embodiment assumes a solid-state, pulse 
counting, segmented detector gamma ray imaging system used in 
a medical imaging application, as shown in Figure 6. The 
detector module 605 has a 4,096 element (channel) segmented 
detector array distributed in a 64 element by 64 element 
rectangular array covering an area of 2 0cm by 2 0cm. Other, 
imaging systems, however, could be alternatively used. 
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[0025] The detector array 605 is subdivided physically into 
32 modules 606, each incorporating 128 contiguous channels. 
Alternatively, modules may be defined electronically rather 
than physically* 

[0026] In the embodiment being discussed, each module such as 
606 is assigned a unique number from 0 to 31. The channel 
numbers in each module are assigned sequentially as 128 
consecutive integers. Channels 0 through 127 are assigned to 
module 0, channels 128 through 255 to module 1, continuing up 
to module 31, to which channels 3,968 through 4,095. In any 
given embodiment, the numbering of modules and channel numbers 
is arbitrary. 

[0027] The signal processor 615 performs an energy analysis 
of the signal from each channel and recognizes as events only 
those pulses corresponding to energy greater than a specified 
minimum energy. It will be seen below that energy filtering 
by this method is useful in identifying DE caused by cosmic 
gamma rays . 

[0028] The techniques disclosed herein are also applicable to 
events discriminated in energy according to other criteria. 
In these embodiments, the logical flow for distinguishing 
between valid data and those associated with distributed 
events may differ from that discussed below. 
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[0029] The gamma camera may employ a serial readout scheme in 
which the channels are read in ascending numerical order, and 
data are output only from those channels that have detected an 
event since the end of the previous readout cycle. Reading a 
channel that has detected a gamma ray requires 2.67|is, and 
reading a channel that has not detected an event requires 
40ns. Therefore, if all channels detect a gamma ray, reading 
out the entire array requires 10.9ms. Other embodiments may 
employ parallel readout of data from the individual channels. 
[0030] In current clinical gamma ray imaging applications, 
the flux of gamma rays emanating from the patient and striking 
the detector array will not exceed 2 00 gamma rays per cm 2 per 
second. At this flux, the average count rate per channel 
corresponds to less than 0.2 events per readout cycle and to 
less than 800 total counts per readout cycle. Distributed 
events such as vibration, electrical discharge, Compton 
scattered cosmic gamma rays and individual gamma rays 
coupling, through electrical parasitics, into neighboring 
channels typically produce false events detected in one 
readout cycle in multiple channels physically located near 
each other, but not necessarily contiguous in channel number, 
or in patterns that are visible in the image but are not 
amenable to being described by one or a plurality of geometric 
relationships. These false events in multiple channels may 
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create a local areal count density over a set of physically 
contiguous channels significantly exceeding the expected count 
density of less than 2 00 events/cm2- second. In addition, 
vibrational DE may produce such elevated local areal count 
densities over more than one successive readout cycle. 

[0031] One occurrence of noise phenomena typically results in 
detection of events in a plurality of channels within the time 
required to read out all of the individual channel signals. 
Therefore events produced by these phenomena are denoted 
distributed events or DE. Distributed events typically cannot 
be discriminated against simply on the basis of primary 
particle energy. Scattering of cosmic rays distributes the 
original energy of the cosmic ray over a plurality of 
interactions with multiple pixels. In the case of distributed 
events caused by vibration or electrical discharge, there is 
no primary particle. There is no specific geometric 
relationship between the locations of individual channels that 
detect individual events caused by a distributed event. 

[0032] Therefore, in SPERIS, it is desirable to discriminate 
against multiple events caused by a DE and not identifiable 
simply by primary particle energy, time coincidence, a 
specific geometric relationship between a plurality of events, 
or by any combination of these criteria. 
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[0033] Distributed events are characterized by time 
concurrence of events in the sense that they cause multiple 
invalid events to be detected during one readout cycle. The 
individual events are not necessarily simultaneous within a 
short time frame of the order of nanoseconds to microseconds. 
In the embodiment described below, discrimination based on 
time concurrence is implicit, because the filter processes 
data acquired during one readout cycle. Other embodiments may 
include explicit discrimination on the basis of time 
concurrence . 

[0034] In many cases, events resulting from distributed 
events cannot be discriminated against on the basis of the 
detected energy of events in individual channels. Although a 
cosmic gamma ray has a well-defined energy, scattering of the 
cosmic ray causes it to release only a portion of its initial 
energy in any single channel. In addition, the cosmic ray may 
not release all of its energy within the detector channels, so 
that summing the energy of all the events associated with the 
cosmic ray typically will not yield the initial energy of the 
ray. The events associated with the cosmic gamma ray, 
however, typically will have an average energy significantly 
greater than the energy of the event of interest or EOI . 
[0035] Techniques are disclosed for detecting false events 
resulting from a distributed event on the basis of the global 
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count density, local count density, average energy of detected 
events, and time period of detection of false events. This 
technique can be implemented in hardware, software or firmware 
within the signal processing subsystem 615 or elsewhere in the 
imaging system. For brevity, this is described as a blink 
filter . 

[0036] The process follows the flowchart, beginning with Fig. 
1. The flowchart of Fig. 1 may be executed in signal 
processor 615. Data read out from a single channel that has 
detected energy that may represent a gamma ray are 
incorporated in a data word. The data word identifies the 
channel number and energy of the event and also includes an 
"ignore" field, e.g., a bit, which can be set to identify 
invalid data words. The data word may also include one or more 
additional fields that can be set to indicate, for example, 
why a data word was treated as invalid. Initially, before 
being processed by the blink filter, each ignore bit is set to 
0 to identify the event as valid. In other embodiments, the 
exact format of the data word, the data incorporated in the 
word, and the valid/invalid identification bit or bits may 
differ from the embodiment. 

[0037] The data words from the individual channels are 
incorporated into a data stream that also includes time words. 
In the embodiment, time words are inserted into the data 
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stream every 0.5ms. In other embodiments of a blink filter, 
one or more of the exact format of the data word, the data 
incorporated in the word, the valid/invalid identification 
code, and other characteristics of the data word may differ 
from the structure of the data word of the embodiment 
disclosed. 

[0038] The data stream from the readout of one detector 
module, including any time words that may be inserted during 
the readout of the module, are assembled into a string. In 
the serial readout of this embodiment, the strings from the 
modules are read out and inserted into the output data string 
in order of increasing module number. The transition from one 
readout cycle, or frame, is identifiable in this embodiment by 
detecting that the channel number of the first data word of a 
later readout cycle will be lower than the last data word of 
an earlier readout cycle. 

[0039] At 100, the raw data stream first is loaded into a raw 
data buffer that holds all the data output from the detector 
module during a time window of W milliseconds. The time 
window of the embodiment is 20ms. In all the figures herein, 
the default values of parameters in the embodiment of the 
example are shown in parentheses. The raw data buffer 
therefore contains data from multiple frames and, as a 
consequence, multiple strings. 
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[0040] Next, a copy of the data in the raw data buffer is 
made at 102 and placed within a copy buffer 105 and all 
further filtering operations employ data in the copy buffer 
created at 105. 

[0041] The data in the copy buffer 105 are first analyzed by 
a raw multiplicity filter at 110. This analyzes the data from 
one module in one readout, denoted a string, based on the 
number, or raw multiplicity, of events in that string. Frames 
in which all strings which have raw multiplicity below a 
threshold value Trm (default value = 2) are assumed to be 
valid, and flow passes to 114, representing no changes being 
made in the corresponding data in the raw data buffer. 100. 
All events in these frames will have an ignore bit of zero 
and, unless the ignore bits of one or more events are changed 
to 1 in subsequent steps, will be included in an image. 

[0042] Frames incorporating one or more strings with a raw 
multiplicity greater than Trm are candidates for analysis as 
distributed events, or blinks, and are further analyzed by a 
global density filter 120. 

[0043] The global density filter 120 determines if the 
average count rate in a frame including a candidate string is 
greater than a high count rate threshold T CR (default value = 
80,000) counts per second. Frames with average count rates 
greater than T CR can be expected to have multiple strings whose 
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raw multiplicity exceeds Trm. Therefore these frames are 
treated as valid, and flow passes to 114, which leaves data 
for these frames are left unchanged in the raw data buffer 
100. Frames with average count rates less than Tar and 
including strings with a raw multiplicity greater than Trm are 
blink candidates and are copied to a candidate buffer 130 that 
holds the data from one frame. 

[0044] Figure 2 illustrates the further analysis of data from 
the candidate buffer 130 using a local density filter shown as 
200. The local density filter 200 determines if the camera is 
in persistence mode, in which data from frames read out over a 
period of T PM (default value = 250) and therefore including a 
plurality of frames/milliseconds are averaged on a channel -by- 
channel basis in the final image. If the camera is in 
persistence mode, images of invalid data from one frame will 
be averaged into the image data for the whole period T PM . 
Therefore, deleting the data from the all the frames stored in 
the Raw Data Buffer 100 will improve the quality of the final 
image. When the camera is in persistence mode, the local 
density filter 200 discriminates against vibrational 
distributed events. 

[0045] Typically, vibrational events will produce a high 
density of false events, and will be associated with a high 
local count rate in one or a plurality of modules. Therefore, 
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they increase the size (number of data words) of the data in 
the candidate buffer 130. 

[0046] When the camera is in persistence mode, the local 
density filter determines if the size of the data in the 
candidate buffer 130 is greater than a threshold of T CB 
(default value = 90) percent of the size of the raw data 
buffer 100 and if the count rate in an individual module, 
denoted the impulse count rate, is greater than a threshold T CR 
(default value =80,000) counts per second. If the camera is 
in persistence mode and both of these criteria are true, the 
data in the entire raw data buffer are marked as invalid at 
202 by setting the ignore bit on all events in the raw data 
buffer. The data in the raw data buffer 100 are transferred 
at 204 to the normal buffer processing system, which excludes 
these data from inclusion in the final image. The raw data 
buffer 100 is then reloaded at 206, and if appropriate, copied 
at 102 into the copy buffer 105. The filter then proceeds to 
analyze the newly loaded data as above. If the data in the 
frame do not meet all of these criteria, the data are then 
analyzed by a mode/buffer size, or MBS, filter 300. In the 
embodiment, rejection of vibrational events in persistence 
mode is performed in all cases. For all other types of 
invalid events, the filter can be turned on or off as desired. 
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[0047] Analysis of data by the MBS filter 300 is shown in 
Fig. 3. If a vibrational DE has occurred, there will be a 
large number of events in the frame and therefore the size of 
the candidate buffer 130 is expected to be large. The MBS 
filter 300 checks to determine if the candidate buffer size is 
greater than a threshold T NP (default value = 90) percent of 
the size of the raw data buffer 100 the data are treated as 
being affected by a vibration. In this case, a string from an 
individual module is likely to contain more invalid events 
than the raw multiplicity threshold Trm. 

[0048] The data in the string are next analyzed to determine 
characteristics that would suggest that the DE was caused by a 
cosmic ray, or was the result of another phenomenon such as 
vibration or an electrical discharge. Distributed events 
caused by all such other phenomena are denoted blinks. 
Because the data in a frame affected by vibration are likely 
to contain multiple invalid events, the threshold T CR for the 
number of events in a string associated with a cosmic ray is 
set to a low value T CM v (default value =4) at 302 to subject a 
high percentage of strings for further analysis. 
[0049] If the data in a frame are unaffected by vibration, a 
cosmic ray striking a module can be expected to cause a 
substantial number of events in a string from that module. 
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Accordingly, the cosmic multiplicity threshold T C r is set to a 
higher value T CM (default value = 10) at 304. 
[0050] All further analysis of data is based on the content 
of strings rather than entire frames. 

[0051] Fig. 4 shows processing the data from the string being 
processed using a blink multiplicity filter 400. Strings with 
a size greater than the threshold T B m (default value = 60) are 
candidates for DE associated with cosmic rays and are assigned 
for processing in an average energy filter 500 shown in the 
flowchart of Figure 5. Strings of shorter length are also 
related to DE caused by cosmic rays and are assigned for 
processing to a cosmic multiplicity filter 410. 
[0052] If a candidate string contains more events than the 
cosmic multiplicity threshold T C r set earlier (see Fig. 3) , the 
string is treated as having been caused by a cosmic ray. The 
string is reported as a cosmic ray at 412 and then processed 
by a filter mode filter 420. If the camera has been set so 
invalid data are to be filtered, the string is marked for 
rejection at 422 by setting to 1 the ignore bits on all events 
in the corresponding string in the raw data buffer. If 
filtering is turned off, the string is treated as valid and no 
change in the raw data buffer at 424. 
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[0053] If the candidate string contains from 0 to T CR events, 
the string is treated as valid and no change is made in the 
raw data buffer at 412. 

[0054] As shown in Fig. 4, candidate strings that have been 
analyzed previously by the blink multiplicity filter 400 and 
have a string size greater than the blink multiplicity 
threshold T B m have been sent to an average energy filter 500, 
whose actions are depicted in Fig. 5. The events in such 
strings have previously been treated as having a cause other 
than vibration. The remaining task of the signal enhancement 
module is to determine if the distributed event was caused by 
a cosmic ray or some other cause, such as electromagnetic 
discharge. Distributed events resulting from such other 
causes are referred to as blinks. 

[0055] Distributed events caused by a cosmic gamma ray 
produce a plurality of individual events with detected 
energies less than the initial energy of the cosmic ray. A 
substantial fraction of such individual events have a detected 
energy significantly greater than the energy of the particles 
being imaged. 

[0056] The average energy filter 500 determines if the 
average detected energy of the events in the string being 
analyzed falls below a cosmic energy threshold of T CE (default 
value = 3 00) keV. If the average energy of the detected 
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events in the string is below T C e# the cause of the events in 
the string is treated as other than a cosmic gamma ray and is 
reported as a blink at 502. Otherwise, the cause of the 
events in the string reported as a cosmic ray. In both cases 
the events in the string are invalid. 

[0057] Following assignment of a cause for the events in each 
string by the average energy filter 500, the string is 
processed by the filter state filter 420. If filtering has 
been turned on, the string is marked for rejection at 422 by 
setting to 1 the ignore bit of all bits in the corresponding 
string in the raw data buffer 100. Otherwise, at 424 no 
change is requested in the raw data buffer 100. 
[0058] When processing of all the data in the raw data buffer 
100 by the signal enhancement module has been completed, the 
contents of the raw data buffer 100 are transferred into the 
normal data processing subsystem, where strings reported as 
being invalid are excluded from the data incorporated into the 
final image. 

[0059] If there are more data to be processed by the signal 
enhancement module, the raw data buffer 100 is refilled and 
processing by the signal enhancement module continues at 102. 
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[0060] A number of embodiments have been described. 
Nevertheless, it will be understood that various modifications 
may be made without departing from the spirit and scope of the 
invention. Accordingly, other embodiments are within the 
scope of the following claims. 

[0061] For example, the above has described the filtering 
being carried out in software in the computer 630, however 
this could also be done in hardware, in firmware, or in other 
processing systems. In addition, while the above has 
described using this technique for gamma radiation, it could 
also be used in any other kinds of energy detection. 

[0062] Also, while the above has described numerous different 
filters, it should be understood that the disclosure 
contemplates different combinations of these filters, which 
may include fewer than all the filters or even a single one of 
the filters being used to set the detection parameter. 

[0063] The embodiment specified one or a plurality of 
criteria for each of the filters. Other embodiments may 
employ none or a subset of the specified criteria and may 
employ one or a plurality of other criteria to perform the 
desired function. In addition, one or more filters with 
multiple criteria may be separated into a plurality of other 

filters. Also, a plurality of filters that collectively 
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perform one or a plurality of functions may be combined to 
form one or a plurality of other filters that will perform the 
same set of desired functions. 
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